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Summary 


The version of FSCBG (Forest Service Cramer Barry Grim) described in this 
technical manual combines and implements mathematical models for aircraft wake effects, 
gaussian line source dispersion, droplet evaporation, canopy penetration, and ground and 
canopy deposition. FSCBG is designed to account for the atmospheric dispersion, 
transport, and deposition of all aerial spray material from the time of release until the spray 
material is either deposited or, in the case of spray drift, until the spray concentration and 
deposition levels become insignificant. Specific calculations made by FSCBG include 
spray concentrations and dosages above and within forest canopies as well as the spray 
deposition within and below forest canopies, resulting from aerial spray releases made 
along single or multiple flight paths. This document serves as the technical resource 
manual for FSCBG. 
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1. Introduction 


The FSCBG (Forest Service Cramer Barry Grim) computer model for predicting 
aircraft spray dispersion and deposition above and within forest canopies results from the 
merging of models for simulating the effect of the aircraft wake on the behavior of released 
spray material, the penetration of spray material into a forest canopy, and determining 
optimum swath widths, application rates, and aircraft altitudes to achieve specific 
objectives. FSCBG combines and implements mathematical models for aircraft wake 
effects, gaussian line source dispersion, droplet evaporation, canopy penetration, and 
ground and canopy deposition. The USDA Forest Service selectively uses aerial spray 
applications to control forest pests while the U. S. Army is interested in assessing the 
effectiveness of chemical and biological defensive strategies. Both agencies are interested 
in achieving a more complete understanding of the behavior of spray material from the ime 
spray is released from the aircraft until it has been deposited or, in the case of spray drift, 
diffused to concentration/dosage levels that are environmentally insignificant. Because 
mathematical spray dispersion models are useful in determining interactions of the many 
factors affecting spray operations, the USDA Forest Service and U. S. Army have 
supported the application and development of these models over the last twenty years. 


Each principal model section of FSCBG is reviewed in subsequent chapters: 
meteorology (chapter 2); evaporation (chapter 3); canopy penetration (chapter 4); near wake 
particle trajectories (chapter 5); and dispersion (chapter 6). This manual is an extension of 
portions of Dumbauld, Bjorklund and Saterlie 1980 and Bjorklund, Bowman and Dodd 
1988. 


2. Meteorology 


All meteorological data are entered in tables in FSCBG, and then curve fits or 
averages are taken to arrive at typical data levels in open terrain or within a canopy. Wind 
speed fits a functional three-layer dependence, while all other variables (temperature, 
relative humidity and wind direction) fit a height-weighted mean value over the depth of 
interest. 


The functional form of the wind speed variation with height depends on the regime 
involved. In open terrain where there is no canopy, or at heights greater than twice the 
height of the canopy, the familiar power law is used. The power-law profile gives the 
wind speed at any height between the surface and the height of the surface mixing layer 


Hm as 
u=u(2)P (2-1) 


where z, is the reference height where the reference wind speed u, is measured, p is the 


power law exponent, z is the height of interest, and u is the resultant wind speed. In 
FSCBG z, is set equal to 10 meters. 


Cionco (1965) developed a wind profile law for describing the variation of wind 


speed within vegetative canopies 


w= uc exp] Oc (7-1 | ; (2-2) 


for any height z upto Hc, the height of the top of the canopy. The parameter Q, is 
called the canopy index, and u, is the wind speed at the top of the canopy. 


A transition region exists (Bergen 1971) between the top of the canopy and twice 
the height of the canopy where the variation of wind speed with height may be given by the 
well-known logarithmic profile 


y= n(2-S*70) (2-3) 


where u’ is the friction velocity, d is the displacement height, and z, is the roughness 
length. In FSCBG d and Z, are set equal to 0.7 H, and 0.14 H, respectively (from 
Dumbauld and Bowers 1983). 


After manipulation, these three equations may be expressed in the form 


f(u) =a + b g(z) (2-4) 


where f(u) is a suitable function of wind speed and g(z) is a suitable function of height. 
In this way the coefficients a and b in Equation (2-4) may be linear least squares fit to 
user-entered data. In FSCBG the user-entered wind speed data in open terrain are fitted to 
Equation (2-1) with the manipulation 


Inu=Inu; +pin( 2) (2-5) 


with linear least squares determining u, and p. The wind speed data within a canopy are 
fitted to Equation (2-2) with the manipulation 


Inu =In Uc - Oe + Oc( FF} (2-6) 


with linear least squares determining ug and G@,. Continuity of wind speed at the top of 
the canopy and at the top of the transition region results in unique coefficients for Equation 
(2-3) with the manipulation 


got jn [Ze )+a° n( 2-870 (2-7) 


and consistency at z=H, and z=2H,. In this way an analytic description of the wind 
speed variation with height is developed from user-entered data. 


If the user enters wind speed data at only a single height, or defines the layer 
average wind speed, Equations (2-1) and (2-2) must be integrated in the vertical direction to 
compute the average wind speed and deduce u, and p in Equation (2-1), and up and O& 
in Equation (2-2). With these values in hand, continuity of the wind speed at z=He and 


z=2H¢ may be invoked to determine u” /k in Equation (2-3). 


Integration of Equation (2-1) for the average wind speed u in open terrain yields 


_U(Hp - 2p )zP(pt+1) 
up = ——_—_—_——_—_——_ 


iru : als (2-8) 


where zp is assumed equal to 2 meters. The value of uy (or u ) then depends on the 
value of the power law exponent p. 


Integration of Equation (2-2) for the average wind speed u within a canopy yields 


DANA MRO Rt 2-9 
ne (1-exp(-a@)) ge 


where O, is determined (Dumbauld and Bowers 1983) by requiring wind speed and its 
derivative to be continuous at the top of the canopy into the transition region (at z= He ). 
This requirement yields 


Qe = bath ina SM dt: = 1.985 


(Hg - 4 r29)1n(He=d 0) 210 


for the values of d and z, assumed in FSCBG. Equation (2-9) then shows that the wind 
speed ug at the top of the canopy equals 2.3 times the mean wind speed within the 
canopy. 


The power law exponent p is determined from tables in Dumbauld and Bowers 
1983 fit by a hyperbolic multiquadric approach (Hardy 1971) to determine smooth first and 
second derivatives with respect to wind speed and net radiation index. Evaluation through 
the table is performed by Gregory square (Barnhill 1977 and Gregory 1974). Bjorklund, 
Bowman and Dumbauld 1984 extend this technique to include evaluation of typical wind 
elevation angle, wind azimuth angle and surface mixing layer height to recover these 
needed meteorological values smoothly. 


3. Evaporation 


In FSCBG the time rate of change of drop diameter due to evaporation is based on 
the expression (Fuchs 1959; Beard and Pruppacher 1971; Goering, Bode and Gebhardt 
1972; and Williamson and Threadgill 1974) 

4 MgDy PA(€s~ &.5) 


+ fl (3-1) 
Mm D pg( PA - es ) M 


Coe 
dt 


where 


Mg = molecular weight of evaporating vapor from the drop (gm/mole) 


= mean molecular weight of the resulting vapor-air mixture, 
approximated by air (Ma gm/mole) 


Dy = molecular diffusivity of the evaporating vapor in air at the drop 
temperature (cm2/sec) 


D = drop diameter (cm) 
PA = airdensity (gm/cm>) 


Pg = density of the drop (gm/cm?) 


€; = partial pressure of the evaporating vapor at the drop surface (mb) 
€oo = partial pressure of the evaporating vapor far from the drop (mb) 
PA = air pressure (mb) 

fy = ventilation factor (no units) 


Air density is calculated from the relationship 


Pa MA (3-2) 


PA= R* Te 


where R”™ is the universal gas constant and Ty is given in terms of the mixing ratio 1,, 
(Beard and Pruppacher 1971) 





T 
rr (1 — | 
t= 1 "0.621974 (3-3) 


+r 
where 
TA =  airtemperature (deg K) 
0.62197 e., 
eS mp acum (3-4) 


The diffusivity D, of vapor into air at the temperature of the surface of the drop Ty; 
depends on the liquid being considered. For water drops FSCBG uses the expression 
(Hall and Pruppacher 1976) 


a ite we) Fo] 3.5 
Dy =0.211| 7 Pa (35) 


where Ty = 273.15 deg K, Pg = 1013.25 mb and Ty is determined from the relation 


LDy Mg{ €s -€,. ] 





Ty=Ta - 3-6 
where (Beard and Pruppacher 1971 and Pruppacher and Rasmussen 1979) 
L == _Iatent heat of vaporization (cal/gm) 
a 
=o) 3 E (3-7) 
i 
a = 0.107 + 0.000367 T; (3-8) 
k = thermal conductivity (cal/cm-sec-deg K) 
1.02 ky \{ “ee 
(snr 2 | 3-9 
kal ki AK (3-9) 


kg = 0.00001 [ 5.69 + 0.017 ( T;- To) ] (3-10) 
ky = 0.00001 [ 3.78 + 0.020(T,-To)] (3-11) 
el ea (3-12) 


The partial pressure eg is defined by the expression (Rasmussen 1978) 


=p 575.0466 + 31.82291 Tj + 1.296028 T;2 (3-13) 
s 93.51611 - Tj 
where 
eo Tea (3-14) 
B = 1 


The vapor pressure e,. is also obtained from Equation (3-13) when Ty is replaced by Ta 


in Equation (3-14) and f is replaced by the relative humidity expressed in fractional form. 
Equations (3-5) through (3-13) are solved by iteration in FSCBG to determine the drop 
temperature, vapor pressure at the drop surface and diffusivity. Then Equation (3-1) is 
solved for the diameter change. 


The ventilation factor fy for water is given (Pruppacher and Rasmussen 1979) by 


fy = 0.78 +0.308 Sc!/3 Rel /2 (3-15) 
for 

1.4<Scl/3 Rel2<514 
and 

fy = 1.0+0.108(Sc1/3 Rel/2)2 (3-16) 
for 


0<Scl/3 Rel/2 << 1.4 
where 


Re = Reynolds number 


VDPA (3-17) 
HA 


Sc = Schmidt number 


eam (3-18) 
Dy PA 


V = _ terminal settling velocity (cm/sec; McDonald 1960) 


absolute viscosity of air (gm/cm-sec; Beard and Pruppacher 1971) 


HA 


0.076342 ( WE (3-19) 
Ta + 296.16 | 296.16 


FSCBG can also theoretically determine the time rate change of drop diameter for 
non-water drops. However, for this option to work, additional input information is 
required. In this case the drop temperature and vapor pressure at the surface of the drop are 
calculated from the expression (Beard and Pruppacher 1971) 


exp(B-C/Tr)-€..=( G5 Arya, |(TA~ TILPA- exp (B- C/Tr)] (3-20) 


where 
Cs = molal concentration of the liquid (moles/cm?) 


and 


B, C = constants obtained from tables expressing variations of vapor 
pressure with temperature (for example, Weast 1981 with some 


manipulation into compatible units) 


Values of k, Dy and L for non-water drops must be specified by the user. When the 


parameters required by Equation (3-20) have been specified, FSCBG uses a Newton- 
Raphson iteration technique to determine the drop temperature, vapor pressure and the time 
rate change of diameter, using the following expression (Fuchs 1959) for the ventilation 
factor 


fy = I Sh (3-21) 


where 


Sh = Sherwood number 


2(1+aSc 13 Re 1/72) (3-22) 


with a =0.3 (Williamson and Threadgill 1974). 


Results obtained using the evaporation model for open terrain calculations are fitted 
to a quadratic by least-squares over the time period required for the drop to evaporate to a 
diameter of 5 micrometers or reach a settling velocity of 0.02 cm/sec , whichever is 
greater. Within canopy calculations are fitted to a quadratic over the time period required 
for the drop to travel from the top of the Canopy to the ground. 


4. Canopy 


The probability that a drop will penetrate a canopy depends on the total number and 
size of vegetative elements encountered (Rafferty and Bowers 1989). If the orientation of 
the vegetative elements is assumed to be random, then the probability of penetration for a 
given path length will be the same for all directions. In FSCBG Pp , the probability of 


penetration, is defined as the probability that a drop traversing a horizontal trajectory will 
penetrate a typical single tree envelope. The value of Pp is assigned for each tree type 


based on foliage density and envelope dimensions. 


In the analytic canopy penetration algorithm the height of the canopy is divided into 
twelve layers and an average tree diameter is computed for each layer. Each tree is 


represented in the kth layer by a cylinder of diameter Wx and height Az, . The average 
length of a drop trajectory through this cylindrical envelope is 


TZ Wy Az (u2 + v2)? 


Lk = eee ee (4-1) 
or Wk V + ux Az | 
where 
Uk = wind speed (horizontal) 
V = _ terminal settling velocity (vertical) 


In the limit of wind speed uy approaching zero, Ly, equals Az, . Similarly, in the limit 
of settling velocity V approaching zero, Ly equals mW, / 4, the mean horizontal path 
length through a circle of diameter W, . A correction factor o may be defined as 


Ly /(% Wi /4 ), designating the number of mean horizontal path lengths traversed by an 
inclined drop trajectory while penetrating a cylindrical tree. Thus, the probability that a 
particle will penetrate one tree is given by 


Py =1-E(1-P} (4-2) 


where 


E = collection efficiency 


10 


The analytic canopy penetration algorithm assumes that the probability of a drop 
encountering a randomly placed tree envelope is equal to the fraction of the horizontal area 
effectively occupied by trees. If the drops fall vertically, this probability is simply the 
fraction of the horizontal area occupied by trees; this fraction must be increased to account 
for horizontal drop transport. The increase in the effective horizontal area of a single tree is 
given by the product of Wk and the downwind distance a drop travels while falling 


through Az, . The resulting probability of encountering a tree envelope is 


Uk Wy Az 
f, = Twi + “k Sk Ak Pt (4-3) 


where py is the tree density (number of stems per unit area). In Equation (4-3) it is 
possible for f, to exceed unity. The overall probability of a drop penetrating a layer is 
then given by 


Pre =[1-Fk(1- Py )] Pik (4-4) 


where Py is from Equation (4-2), and I, and Fy are the integer and fractional parts of 
f, , respectively. For example, if fk equals 2.4, I, equals 2 and Fy equals 4. In 
this case, six-tenths of the drops penetrate two trees and the remaining four-tenths penetrate 
three trees. 

The collection efficiency E is determined by impaction and interception with 
cylindrical representations of the tree elements. Two parameters are evaluated 


St = Stokes number 


2 
mee Dis U (4-5) 


9 HAD; 


9 = 22k DU (4-6) 
HA Pg 


where 


Pg = density of the drop (gm/cm>) 
D = drop diameter (cm) 
U = particle speed (cm/sec) 


11 


HA = absolute viscosity of air (gm/cm-sec) 


D,; = characteristic size of a tree element (cm) 


PA = airdensity (gm/cm?) 


St and @ are then used with the cylindrical impaction curves of Golovin and Putnam 1962 
to determine the impaction efficiency E; (Figure 4-1). This efficiency is augmented by an 


interception effect (Pulley and Walters 1990), which for a cylinder gives 


E=Ri+Ch 


with 


C= 1+ 0.33 St 0-96 


Impaction Efficiency 


Figure 4-1. 
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Cylindrical impaction efficiency curve for @ = 0 (Golovin 


and Putnam 1962). 
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5. Near Wake 


A Lagrangian approach is used to develop the equations of motion of discrete 
particles released from aircraft, with the resulting set of ordinary differential equations 
solved exactly from step to step. The formulation of the near-wake equations, particularly 
the ensemble averaged approach used to develop turbulent properties of the particle- 
atmospheric interaction, is detailed in Bilanin and Teske 1984. The near-wake 
implementation is consistent with AGDISP (Teske 1990). 


Released material is assumed to be spherically shaped. Particle flight path as a 
function of time after release is computed as the locations (X,Y,Z) for all particles 
included in the simulation; velocity is denoted by (U,V,W). The interaction of the 
material with the turbulence in the environment creates turbulent correlation functions for 
the position and velocity, < yv > and <zw > ; for the velocity variance, <vv> and < 
ww > ; and for the position variance, <yy > and <zz>. The square root of these last 
two variables gives the horizontal and vertical standard deviations of the material motion 
about the mean described by Y and Z. 


EQUATIONS OF MOTION 


In the near wake, released material is assumed to be influenced by two forces: 
weight and aerodynamic drag. The governing equations are derived in Bilanin and Teske 
1984 and repeated here for completeness 





Mean Equations: 
dX; 1 
=I 0 raat 5-1 
AD le 1 oa 1 (5-1) 
Si =; (5-2) 
Turbulent Correlations: 
= < XjxXj > = 2 < xjvj > (5-3) 
Te oat ae (5-4) 
a cyivi>=2(<ujvj>-<vjvj,>){L P 
dt ViVi ( ivi ivi Neo (5-5) 


Ds 


In Equations (5-1) to (5-5), Xj, Vj; and Uj; are the ensemble averaged ith 
components of material position, material velocity and local fluid velocity, respectively, 
while x;, v; and u; are the fluctuating ith components of material position, material velocity 


and local fluid velocity, respectively; t is time and g; is the gravity vector. 


Inherent in the equations is a relaxation time which is essentially the e-folding 
time for the released particle to come up to speed with the local fluid velocity (for V; to 
approach and equal Uj ). The relaxation time is defined as 


(5-6) 


where D is the diameter of the drop, py is its density, pa is the density of air, V,,) is 


the relative velocity between the material velocity and the local background velocity, and 
Cp is the particle drag coefficient. The semi-empirical formula for Cp (Langmuir and 


Blodgett 1946) is 


Cp = ae 1 + 0.197 Re9-63 + 0.00026 Re!.38 ) (5-7) 


where Re is the Reynolds number defined as 


Rep relen sy (5-8) 
HA 
with 
Vre] = Telative velocity = IU; - V;| 
HA = absolute viscosity of air 


Equations (5-1) to (5-5) cannot be solved without specifying relationships for the 
quantities <xju; > and <ujv; >, the correlations of the particle position and particle 
velocity fluctuations, respectively, with the local fluid velocity fluctuation. These 
expressions are developed by integrating their ensemble averaged frequency spectra using a 
spectral density function for transverse velocity fluctuations in isotropic turbulence (von 
Karman and Howarth 1938). The results (detailed in Bilanin and Teske 1984) give 
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<xjuj>=2| ey K+ %] (oA, 


<uivi>= gq? K (5-10) 
Lb aa5 a 
with 


q2 = mean square turbulence level 


= <uu>+<w>+<ww> 


p\*|[,.% Elie 
Tr 


eek sn 
Gay 


and 7, is the travel time of the particle through a turbulent eddy of scale A , adjusted for 
the passive tracer limit 


A 


With the position and velocity information available for all released material at any 
time during the simulation, Equatiors (5-1) to (5-5) may be integrated exactly for the 
solution at the next time step (Teske '“88). The mean Equations (5-1) and (5-2) may be 
combined to give 


2x. vu; 
CEES «SS ES Oe (5-13) 


di2 dt 
An exact solution of this equation, from a time of 0 toa time of t, yields 
Xj = Xo - T (Uj - Vo + Bi tp) 1-exp(-t/t) + (Uj+g) %)t (5-14) 


Vj =- (Uj - Vo + 8) Mp ) exp (- t/t) + (Uj + gj tT) (5-15) 
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where 


Xo = initial particle position 


Vo = initial particle velocity 


Similar equations may be written and solved for the turbulent correlations. Any step result 
Serves as initial conditions for the next step integration, etc., through the entire simulation. 
The assumption that the background conditions Uj, <xjuj > and <ujvj > are constant 


across a time step permits solution of Equations (5-1) to (5-5) unrestricted by the size of the 
released material. Computational solution times are reduced significantly, and any drop 
size may be tracked by the equations. 


FLOW FIELD MODELING 


The behavior of the released material is intimately connected to the local background 
mean velocity U; and turbulence field q? through which the material is transported. The 


near wake module in FSCBG contains a number of simplifications for the flow field 
velocity and turbulence levels additionally imposed in the aircraft wake. These models are 
summarized below. 


Fixed-Wing Rolled Up Tip Vortices 


When an aircraft flies at a constant altitude, the aerodynamic lift generated by the 
lifting surfaces of the aircraft equals the aircraft wei ght. The majority of the lift is carried 
by the wings, and generates one or more pairs of swirling masses of air (called vortices) 
downstream of the aircraft. If the roll-up of this trailing vorticity can be approximated as 
occurring immediately downstream of the wing, then the mean velocity field that results 


may be simply characterized by the aircraft semispan s , circulation TF and load 
distribution. For a rectangularly loaded wing 


ie 
2s pau (5-16) 
where 


W, = aircraft weight 


aircraft flight speed 


= 
8 
Tl 
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Vortex Circulation Decay 


The flow field generated by a simply loaded wing includes four vortices: the two 
from either wing tip, and image vortices below the surface. The local ambient velocity 
from the vortices is the vector sum of the four swirl velocity components, each written as 


Vs = 5h (5-17) 


with its direction perpendicular to a line connecting each vortex centroid and the 
observation point (the distance r). The observation point might be a particle position or 
the location of another vortex (which also moves because of the other vortices present in 


the flow). An inviscid vortex of constant positive circulation strength I\, descends toward 
the surface while increasing the separation distance between its companion vortex. 


In the atmosphere, turbulence acts to decay the vortex strength. A simple decay 
model, developed in Donaldson and Bilanin 1975, results in a vortex decay of the form 


r=T>% exp { - Bat] (5-18) 


For vortices out of ground effect, a typical value for B = 0.41 : while a detailed 
examination of Program WIND data in ground effect (Bilanin et al. 1989) gives 
Bq = 0.56 m/sec. 


Helicopter In Forward Flight 


The helicopter model includes both hover downwash and the tip vortex pair by 
partitioning the helicopter weight between the two effects as a function of time. The hover 
downwash model is taken from actuator disk theory for a propeller and may be written as 


F W;=2 pa R2w? (5-19) 
while the strength of the vortex pair may be found from 


(1-F)W,=2paRU,,T (5-20) 


where 


17 


W, = helicopter weight 


Wq = downwash at the rotor plane 
R= = radius of the rotor 
[I = tip circulation strength 


The function F is determined from the expression 


Feexp|-kATot (5-21) 
onR2 
where 

Ht = forward advance ratio 

© = solidity 


Ig = solution to Equation (5-20) for F=0 


The constant k relates helicopter roll-up (around a circumference at the blade tip) to fixed- 
wing roll-up. Predicted results from a helicopter wake model (Bliss, Teske and 
Quackenbush 1987) suggest that the downwash flow field quickly transitions (within two 
blade diameters) into a vortex pair. A typical value for k=7. 


The tip vortices influence the flow field everywhere, while the downwash affects 
only the spreading region beneath the rotor and the fluid column directly above the rotor. 
Within the boundaries of the rotor blade and the dividing streamline, the velocity is 
assumed to follow a stagnation point flow, while the vortex pair created by the helicopter is 
assumed to follow the dividing streamline. 


Material released ahead of the helicopter (spray boom forward) is assumed to 
encounter a streamline pattern similar to flow around a circular cylinder (the actuator disk 
downwash, Kuethe and Schetzer 1964) 





Zin 
ib 1-Bo4i0y (5-22) 
2U,,R2xy 
oceania (5-23) 


18 


until the material crosses the plane of the helicopter shaft centerline. 
Jet Engine 


_ The jet engine exhaust is modeled as a turbulent circular jet using the similarity 
analysis reported in Schlicting 1968. The results give the axial and radial velocities as 


avon 


Naxial ~ 9 ex Tare (5-24) 
(8) 


4 


AE (v2 


ion (5-25) 


2 \2 
i 
(142 


vata 4 fF 
with 


Alan 
n=l fa BAG (5-26) 


€=0.0161 ,/ me (5-27) 


where T is the jet thrust and x is the downstream distance corrected for the virtual origin. 
The radial velocity v,,4j,] becomes part of the local fluid velocity, while the axial velocity 


Uaxial 1S also used to evaluate the turbulence level within the jet exhaust 





q2 = 0.2034 u2 (5-28) 


where the proportionality constant is determined from the centerline decay of turbulence in 
a free jet (Wygnanski and Fiedler 1969). 


L9 


Propeller 


The propeller is modeled as an actuator disk, where the incremental velocity AU 
above the flight speed U,, is related to the actual thrust produced by the propeller 


T=2pa R2 AU(U,, +AU | ? (5-29) 


where R is the radius of the propeller. In steady flight the thrust equals the aircraft drag, 
so that 


T=Cp1p, U2 s (5-30) 


where S is the wing planform area and Cp is the aircraft drag coefficient. Combining 
Equations (5-29) and (5-30) to eliminate thrust, results in 


AU =1{-14 4/1428 | 
Pes 2 mR2 (5-31) 


An alternate expression for €, defined in Equation (5-27), may be obtained from 
Schlicting 1968 as 


€ = 0.0256 R AU (5-32) 


Combining this equation with Equation (5-27) gives the effective thrust level for the 
propeller as 


Teff = 2.5275 pa R2 AU2 (5-33) 


to permit the use of Equations (5-24) to (5-28) for the propeller when Equation (5-33) is 
substituted for T. The swirl velocity generated by the propeller is assumed to be linear in 
R_ out to its maximum value, and then zero for larger values of r. The resulting 
integration for swirl energy yields 
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us cps 


wn. OR3(U,, +AU] ae 


Vtip 


where 1) is the propeller efficiency and Q is the propeller rotation rate. 
DEPOSITION MODELING 


As material approaches the surface, deposition begins, and continues until all 
material is deposited (if evaporation occurs, some of the material will be left in the 
atmosphere to drift). Ground deposition is computed by assuming that the concentration of 
material around the mean may be taken as gaussian 





24 ps 
ies #4] | -@-¥) | peizeZ) | 
5a exp aD exp wee (5-35) 


where the released material is at position (Y,Z) , and the standard deviation is o. Thus, 
as the material approaches the surface, it deposits incrementally, and continues depositing 
at its approach rate after impact. Equation (5-35) is integrated from far below the surface to 
the material location to give 


olay ea 4") Ze | 
Mplane aks aa erfo( A} (5-36) 


where erfc is the complementary error function. Ground deposition is obtained by 
summing «.'' incremental contributions to Mplane as the integration proceeds. It may be 
seen that io; material falling vertically toward the surface, the ground deposition pattern 
generated by Equation (5-36) will be identical to a gaussian deposition pattern. The use of 
Equation (5-36) leads to more realistic ground deposition distributions when the particles 
are traveling with nonzero horizontal velocities. 
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6. Dispersion 


The spray dispersion models incorporated into FSCBG are similar to aircraft spray 
models described by Dumbauld, Rafferty and Cramer 1976. These models were developed 
previously for the U. S. Army Dugway Proving Ground and are based on the generalized 
diffusion model concepts described by Cramer et al. 1972. These models were also used 
in the CBG computer program, with their mathematical basis given in Appendix B of the 
report by Dumbauld, Rafferty and Bjorklund 1977. These models calculate peak 
concentration, dosage and deposition downwind from nearly-instantaneous elevated line 
sources oriented at an arbitrary angle with respect to the mean wind direction. The axis of 
the spray cloud is assumed to be inclined from the horizontal plane by an angle that is 
proportional to Vj /u , where Vj is the gravitational settling velocity for the jth drop size 
Category and wu is the mean transport wind speed (above the canopy). When evaporation 
is negligible, this inclination angle is invariant with distance from the line source. For 
evaporating drops the angle changes with distance from the source because V; depends on 


the size of the drop. It is assumed that drops dispersed upwards by turbulence are totally 
reflected downwards at the top of the surface mixing layer, but the fraction of drops 
reflected at the ground is a variable input parameter for each drop size category. The 
models use a Cartesian coordinate system for a line source of length L at a release height 


H and a calculation point at R(e,5,z) as shown in Figure 6-1. 
VOLUME SOURCE MODELS 


With these assumptions, the dosage for an instantaneous volume source is given by 
the expression 
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and 


where 


Hm 


(6-2) 





Strength of the volume source 
standard deviation of the crosswind spray distribution 


Standard deviation of the vertical spray distribution 


the alongwind, crosswind and vertical coordinates of the point at which 
the dosage is calculated 


mass fraction of the total source strength in the jth dropsize Category (total 
number of dropsizes is NJ ) 


reflection coefficient for the median drop by mass in the jth dropsize 
category (Figure 6-2) 


depth of the surface mixing layer beneath a capping inversion 


The lateral and vertical growth of the spray cloud due to turbulent mixing is 
assumed to be rectilinear. The expressions for 6, and oO, then become 


y 


Oy = OA (X+Xy) (6-3) 


OA 


OF 


OA (X+Xy) 


O7 = Of (xX+xy)= k (6-4) 
= standard deviation of the wind azimuth angle (in radians) 
= standard deviation of the wind elevation angle (in radians) 
Oo 
= 4 (6-5) 
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as 
| 


k O, 
virtunt distances = -xp =20-xp (6-6) 


Mf OF 
and 
Og = Standard deviation of the cloud distribution at the distance xp 
XR = distance downwind from the volume point source 


The amount of spray material deposited on the ground or entering the forest canopy 
through gravitational settling is obtained from the expression 


NJ 0 
= eV asa e lem Oy Ae ndely Sal naet| Wana (sey ag i 
2 V20 Gy aig 23, x bias ¥20 Oz J. bee Ce 


After substitutions for o, and O, into Equation (6-7), and performing the indicated 
integration and differentiation, the deposition equation takes the form 
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LINE SOURCE MODELS 


The expressions for dosage and deposition downwind from a line source oriented at 
an arbitrary angle 0 with the wind direction are derived through consideration of the line 
source geometry shown in Figure 6-1. A finite line source of length L is directed along 
the 6 coordinate at height H with one end of the line source at the point e=0, 6=0 
and z=H. From the geometry shown in Figure 6-1, it may be seen that 





x =x'- 6' sin® (6-11) 

MESICOSTE © tan OF 5 (6-12) 
cos 8 

x'=e€cos8+6 sin @ (6-13) 


When Equations (6-11), (6-12) and (6-13) are substituted for x , y and x’ in Equations 
(6-1) and (6-8), the dosage and deposition at any point R downwind from the line source 
becomes 


lL 

DOS] = | DOSy dé’ (6-14) 
0 
p 

DEP; = | DEPy dé' (6-15) 
0 


The solution procedure is laborious, but leads to the following expressions. 
DOSAGE ABOVE CANOPY 


The dosage (mass times time divided by volume) above the canopy is the sum of 
contributions from drops in the dispersing cloud and from vapor produced by drop 
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evaporation. The above canopy dosage for the contribution from drops is given by the 
expression 
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where 
a = Y2oa(x'+x,} (6-17) 
b= Yoa(x'+xy-Rsine) (6-18) 
V: 
CMe yu iy (6-19) 
V: 
D = 2iHm+H+ 27 (6-20) 
N = (x' +x, )cot 6 +x’ tan @ - 5 (6-21) 
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(6-22) 


(6-23) 


(6-24) 


(6-25) 


(6-26) 


(6-27) 


(6-28) 


(6-29) 


(6-30) 


(6-31) 


(6-32) 


As inspection of Equation (6-16) and the definitions of the parameters appearing in 
the equation show that the model cannot be used when the wind direction is exactly parallel 


to the line source ( @ = 90 degrees) or when the wind direction is exactly perpendicular to 
the line source ( 6 = 0 degrees). For the case when 0 = 90 degrees, FSCBG uses the 


exact solution by making the following observations (from Figure 6-1 with rotation of the 
wind to along the line source) 


X= Xx: - 6 
yee 
x'=6 
(6-33) 
cos 8 =0 
N=e 
sin 8 = 1 


and defining x' as the downwind distance from the upwind edge of the line source and y 
as the crosswind distance from the line source (perpendicular). Equation (6-16) then 
applies with these substitutions. 


For the case @ = 0 degrees, FSCBG also uses the exact solution, with 


X=c 
y=6'-8 (6-34) 
X =e 


where x is now the downwind (perpendicular) distance from the line source and y is the 
Crosswind distance from the center of the line source. In this limit the standard deviations 


of the crosswind distribution Oy and vertical distribution OG, are defined by the 
expressions 


Oy =OAX+05 (6-35) 


Oz = OE X+ Oo (6-36) 
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The result becomes 


NJ 
_-_ QW AE LW2+y)\, | L2-y 
DOS, Be ont =, 93|%| 2 oy He Sy A(x,z) (6-37) 








Vapor emitted from evaporating drops also contributes to the total dosage above the 


canopy. In FSCBG the evaporation model determines the mass of vapor emitted in Ah 
height units along the cloud axis as it descends to the canopy. The distance downwind 
from the line source where the vapor is emitted is given by ut where t is the time from 
the evaporation model when the cloud axis passes through the midpoint of the height 


interval Ah. The source dimension of the vapor cloud emitted over the Ah height interval 
is defined by 


See (6-38) 


under the assumption that the vapor is rectangularly distributed at the point of emission. 
The vapor dosage contribution is then calculated from the appropriate equations with 


Vj = 0, Ns 1, and xp =ut. The total dosage at each calculation point is determined by 
summing the contributions from drops and vapor. 


DEPOSITION ABOVE AND WITHIN CANOPY 


The deposition (mass divided by area) above the canopy is given by the expression 


NJ exp| G2.-p 
pen = 28 § (1-y) olf | 
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with the same substitutions as before, in addition to 


Vix 
B=H+—— 


E=k2D2+N2 
F=k2 B2+N2 
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(6-39) 


(6-40) 


(6-41) 


(6-42) 


(6-43) 


(6-44) 


(6-45) 


(6-46) 


The deposition model that applies when the wind direction is exactly parallel to the 
line source may be obtained from Equation (6-39) using the substitutions defined by 
Equation (6-33). Deposition for the case when the wind is exactly perpendicular to the line 


source is found using the substitutions defined by Equation (6-34), resulting in 
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Again, x is the downwind distance from the source and y is the crosswind distance from 
the center of the line source; Oy and O, are defined by Equations (6-35) and (6-36). 


DOSAGE WITHIN CANOPY 


In calculating within-canopy deposition, FSCBG makes the simplifying assumption 
that the effects of within-canopy diffusion are negligible. Thus, drops that penetrate 
through the canopy are assumed to follow simple ballistic trajectories. Deposition within 
the canopy is related to deposition above the canopy by a canopy penetration algorithm that 
determines the fraction of drops in each size category that reach a given height within the 
canopy (Section 4). 


To compute dosage within a canopy, Bowers 1991 suggests that with simplifying 
as:umptions 


D 
DOS) = ii (6-48) 


EVAPORATION EFFECTS WITHIN CANOPY 


When no evaporation occurs, the calculation of deposition using Equation (6-39) is 
relatively straightforward because the gravitational settling velocity V; is invariant with 
time. However, when evaporation occurs, both dropsize and gravitational velocity vary 
with travel distance from the source, and the deposition calculations become more 
complicated. Within a canopy FSCBG follows a stepwise procedure for each jth dropsize 
as summarized below. 
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1. FSCBG uses the above-canopy (open terrain) evaporation model described in 
Section 3 to establish a quadratic curve relating the gravitational settling velocity Vj of the 


median drop in the jth dropsize category as a function of travel distance x from the source 
for above-canopy meteorology. 


2. The above-canopy evaporation model is also used to calculate the height 
trajectory of the jth dropsize as a quadratic function of travel distance. 


3. The within-canopy evaporation model is used to calculate the trajectory of the 
median drop in the jth dropsize category, using within-canopy meteorology and assuming 
that the drop entered the canopy with the settling velocity Vie given by the above-canopy 


evaporation model at the canopy height. In addition to calculating the mass loss within the 
canopy, the within-canopy evaporation model calculates the travel distance of the median 
drop from the time it enters the canopy until it strikes the ground. 


4. FSCBG internally develops a curve relating the difference between the 
receptor distance from the source and the travel distance in the canopy for all dropsize 
categories versus the settling velocity at the canopy top Vic for each jth category. 


5. The intersection of the curve developed in Step 4 and the above-canopy 
evaporation curve for the jth dropsize category defines the distance Xef¢ from the source to 


the point where the jth median drop would have to enter the canopy in order to reach the 
calculation point. 


6. The effective height Herp of the cloud centroid for the jth dropsize category 


at the distance from the source determined in Step 5 is obtained from the height versus 
distance curve determined in Step 2. 


7. _ The deposition model equations yield the corrected deposition values when an 
effective value of Vj is used. This effective value is found from the expression 


Pe 


The calculations then proceed in the usual fashion except for appropriate adjustments for 
mass loss due to canopy impaction. 


PEAK CONCENTRATION 


The peak concentration (mass divided by volume) which occurs at the calculation 
point is determined from the expression 


DOS] 1 
CON|p = ae (6-50) 
Xx 
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where ©, is the standard deviation of the cloud in the alon gwind direction defined by 


ox -| ou ce lk (6-51) 
and 
L{x) = alongwind cloud length 
: O6Au(,’.Lsing | (6-52) 


Au vertical wind-speed shear in the layer containing the cloud 


Peak concentrations obtained using FSCBG are strictly correct only for crosswind releases. 
Errors made for alongwind releases and releases at an angle to the wind direction decrease 
at distances beyond the more downwind end of the line source. 






Wind 
Direction 





Figure 6-1. Line source geometry. 
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Settling Velocity (m/sec) 





“0.0 0.2 Ome G 0.8 1.0 


Reflection Coefficient 


Figure 6-2. Relationship between the gravitational settling velocity Vj 
and the reflection coefficient Y; at the ground (from 
Dumbauld, Rafferty and Cramer 1976). 
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